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The  water  motion  in  a  simplified  serpentine  channel  is  investigated  to  study  the  dynamic  water  behavior 
in  the  serpentine  channel  in  a  proton  exchange  membrane  fuel  cell  (PEMFC).  The  effects  of  corner  shapes, 
the  number  of  water  inlet  pores,  the  air  inlet  velocity  and  the  water  inlet  velocity  on  the  removal  of  water 
are  discussed.  Results  show  that  the  water  drains  out  more  quickly  in  a  channel  with  a  round  corner  than 
in  a  channel  with  a  sharp  corner.  A  channel  with  hydrophobic  walls  has  both  a  higher  water  removal 
frequency  and  higher  water  coverage  on  the  bottom  wall  than  does  a  channel  with  hydrophilic  walls. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Water  management  is  of  great  importance  for  a  PEMFC  oper¬ 
ating  at  the  normal  temperature  (80  °C).  Low  water  content  in  the 
membrane  causes  low  proton  conductivity  and  consequent  the  bad 
performance  of  a  PEMFC,  while  too  much  water  leads  to  “flood¬ 
ing”  that  will  block  the  diffusion  of  the  reactant  gas  and  eventually 
worsens  the  performance  of  a  PEMFC.  Now  more  and  more  research 
works  focus  on  the  dynamic  behavior  of  water  in  a  PEMFC,  espe¬ 
cially  water  motions  in  channels. 

There  is  an  abundant  body  of  literatures  on  water  management 
of  PEMFCs  in  both  experiments  [1]  and  numerical  simulations. 
Researchers  utilize  numbers  of  techniques  to  investigate  the  water 
transport  inside  fuel  cells,  including  direct  optic  observation  [2-9], 
neutron  imaging  [10-13],  X-ray  radiography  [14,15],  MRI  [16],  THz 
imaging  [17],  02-sensitive  dye  films  [18]  and  water  sensitive  paper 
[19].  These  researches  mainly  focus  on  the  process  in  which  water 
generates  from  the  gas  diffusion  lay  (GDL)  to  channels  and  in  which 
water  moves  inside  cathode  channels  and  provide  a  good  under¬ 
standing  for  the  water  behavior  inside  fuel  cells. 

Numbers  of  modeling  works  mainly  focus  on  the  water 
management  of  fuel  cells  [20-27]  and  the  GDL  [28-35].  Some 
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other  modeling  works  investigated  the  motion  of  given  water 
droplets/films  in  channels.  Quan  et  al.  studied  different  given 
water  droplet/film  distributions  under  different  fuel  cell  operat¬ 
ing  conditions  to  obtain  the  water  behavior  inside  a  serpentine 
micro-channel  [36].  Five  cases  with  different  initial  water  distri¬ 
butions  in  a  serpentine  channel  are  studied.  Zhan  et  al.  studied 
the  effect  of  hydrophilicity  of  a  graphite  plate  and  the  hydropho¬ 
bicity  of  a  GDL  on  the  motion  of  a  given  water  droplet/film  in  a 
straight  channel  and  a  serpentine  channel  [37].  Results  show  that 
the  hydrophilicity  of  channel  and  the  hydrophobicity  of  the  GDL 
play  an  important  role  in  the  liquid  water  motion.  Jiao  et  al.  inves¬ 
tigated  air-water  flow  behaviors  inside  serpentine  flow  channels 
[38].  Cases  with  five  spherical  water  droplets  and  water  films  with 
a  thickness  of  0.2  mm  attached  to  different  places  in  channels  are 
studied.  Result  shows  that  water  distribution  and  pressure  drop 
vary  in  these  five  cases.  Akhtar  and  Kerkhof  studied  the  dynamic 
behavior  of  liquid  water  transport  in  the  cathode  with  a  tapered 
channel  design  [39].  The  motion  of  a  liquid  water  droplet  of  spher¬ 
ical  shape  with  a  volume  of  0.268  mm3  suspended  in  the  center 
of  the  channel  is  investigated.  Results  show  that  hydrophobic  top 
and  side  channel  walls  and  a  hydrophilic  bottom  wall  are  helpful 
to  remove  water.  Kim  et  al.  studied  the  effects  of  geometries  and 
surface  properties  of  flow  channels  on  the  liquid  water  exhaust 
capabilities  [40].  Motions  of  two  water  droplets  with  a  total  vol¬ 
ume  of  1 .0  |ul1  assumed  to  initially  exist  inside  the  flow  channels  are 
studied.  Results  show  that  hydrophobic  channels  exhausted  liquid 
water  much  faster  than  did  hydrophilic  channels,  but  the  stagnant 
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liquid  volume  in  corners  of  hydrophobic  channels  was  larger  than 
that  of  hydrophilic  channels.  Le  and  Zhou  investigate  the  behav¬ 
ior  of  water  droplets  in  serpentine  channels  [41  ].  Motions  of  eight 
water  droplets  initially  attached  in  the  inlet  part  of  channels  are 
studied.  Results  show  that  the  flow  is  unevenly  distributed  and  that 
small  droplets  remain  in  the  corners  of  the  channel.  Le  et  al.  also 
indicate  that  the  eight  water  droplets  affect  the  motion  of  each 
other.  Above  works  made  an  insight  in  the  water  motion  in  chan¬ 
nels  of  a  PEMFC.  However,  these  works  focus  only  on  the  behavior  of 
pre-existed  water  droplets/films  in  channels,  without  considering 
the  consecutive  emergence  of  water  from  the  GDL  to  channels. 

Recently,  studies  on  the  motion  of  water  generating  from  the 
GDL  are  carried  out.  Quan  and  Lai  studied  effects  of  the  hydrophilic- 
ity  of  channel  surfaces,  the  channel  geometry  and  the  air  inlet 
velocity  on  water  behaviors  in  a  channel  in  which  water  gener¬ 
ates  from  the  whole  surface  of  the  GDL  [42].  Results  show  that 
compared  with  a  round  corner  channel,  a  sharp  corner  channel 
could  be  a  better  design.  Ding  et  al.  investigated  gas-liquid  two- 
phase  flow  patterns  in  cathode  channels  [43].  The  water  motion 
in  a  channel  with  320  pores  on  the  GDL  surface  with  the  same 
diameter  of  400  [xm  is  simulated.  Results  show  that  using  a  more 
hydrophobic  GDL  surface  is  helpful  to  expel  water  from  the  GDL 
surface,  but  may  also  increase  the  pressure  drop.  Zhu  et  al.  studied 
the  dynamic  behavior  of  liquid  water  entering  a  straight  channel 
through  a  GDL  pore  via  a  2D  model.  Effects  of  channel  size  and  pore 
size  are  taken  into  account  [44,45].  On  second  thoughts,  Zhu  et  al. 
studied  the  dynamic  behavior  of  liquid  water  entering  a  straight 
channel  through  a  GDL  pore  via  a  3D  model.  Effects  of  static  con¬ 
tact  angle,  air  flow  velocity,  water  injection  velocity  and  dimensions 
of  the  pore  are  studied  [46].  Results  show  that  hydrophobic  walls 
result  in  earlier  detachment  of  water  droplets,  while  hydrophilic 
walls  results  in  spreading  of  the  water  droplet  on  the  GDL  sur¬ 
face.  Le  et  al.  studied  liquid  water  behaviors  in  a  proton  exchange 
membrane  fuel  cell  cathode  with  serpentine  channels  [47].  Results 
show  that  key  factors  that  determined  the  shapes  and  behaviors 
of  the  droplets  are  the  velocity  of  the  air  flow,  the  contact  angle, 
the  surface  tension,  and  the  viscosity  of  the  droplet.  Ding  et  al.  also 
studied  water  motions  in  five  different  water  injection  structures 
[48].  Results  show  that  more  hydrophobic  GDL  surface  and/or  more 
hydrophilic  channel  walls  would  be  helpful  to  remove  the  liquid 
droplets  from  the  GDL  surface  to  the  channel  walls.  These  works 
focus  on  the  behavior  of  the  water  emerging  from  the  GDL  and 
provide  a  good  understanding  of  the  dynamic  water  behavior.  How¬ 
ever,  these  works  focus  on  the  water  behavior  in  a  straight  cathode 
channel  with  only  one  water  injection  pore  or  water  behavior  in  a 
channel  with  a  whole  surface  generating  water. 

In  the  present  work,  we  aim  to  investigate  the  effect  of  water 
emerging  from  one  pore  on  the  formation  of  and  subsequent  the 
behavior  of  water  emerging  from  other  pores  and  the  effects  of  the 
hydrophobicity/hydrophilicity  of  channel  walls,  the  shape  of  the 
corner,  the  air  inlet  velocity  and  the  water  inlet  velocity  on  the 
water  behavior  and  to  investigated  the  volume  fraction  of  water 
inside  the  channel  chamber  and  the  facet  fraction  of  water  on  the 
bottom  wall  to  study  the  mechanism  of  the  water  transport  in  ser¬ 
pentine  cathode  channels. 


2.  Numerical  models 

2.1.  VOF  model 
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Fig.  1.  Schematic  of  the  computational  domain. 


fraction  of  one  (or  more)  of  the  phases.  For  the  gth  phase,  this 
equation  has  the  following  form  [49]: 
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where  a ,  p,  v,  Saq,  mpq  and  mqp  are  volume  fraction,  density,  veloc¬ 
ity,  source  term,  mass  transfer  from  phase  q  to  phase  p  and  mass 
transfer  from  phase  p  to  phase  q,  respectively. 

The  momentum  equation  is  shown  below: 

+  V  .  (pw)  =  - Vp  +  V  ■  [/x(Vv  +  VDr)]  +pg  +  F  (2) 

where  /z,  p,  F  and  T  are  stand  for  dynamic  viscosity,  pressure,  source 
term  and  temperature. 

The  energy  equation,  also  shared  among  the  phases,  is  shown 
below: 


|(p£)  +  V  .  (5(pE  +  p))  =  V  .  (keffVT)  +  Sh  (3) 

where  /<eff,  Sh  and  E  =  ^q=1aqpqFq/^g=1a'qpq  are  effective  ther¬ 
mal  conductivity,  source  term  and  energy. 

Surface  tension  is  a  force,  acting  only  at  the  surface,  that  is 
required  to  maintain  equilibrium  in  such  instances  and  is  shown 
below: 

P2-Pi=a(T  +  T)  (4) 

where  pi  and  p2  are  the  pressures  in  the  two  fluids  on  either  side 
of  the  interface,  a  is  the  surface  tension  coefficient,  Ri  and  R2  are 
the  radii  in  orthogonal  directions. 

The  force  at  the  surface  can  be  expressed  as  a  volume  force  and 
it  is  the  source  term,  which  is  added  to  the  momentum  equation.  It 
has  the  following  form: 


fVol 


p/qVcq 
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where  k  =  V  •  (n/|n|)  is  surface  curvature. 


(5) 


The  VOF  model  is  used  to  model  liquid  and  gas  fluids  by  solving  a 
single  set  of  momentum  equations  and  tracking  the  volume  fraction 
of  each  of  the  fluids  throughout  the  domain. 

The  tracking  of  the  interface(s)  between  the  phases  is  accom¬ 
plished  by  the  solution  of  a  continuity  equation  for  the  volume 


2.2.  Geometric  model 

Fig.  1  shows  a  schematic  illustration  of  the  computational 
domain  and  corresponding  mesh.  In  the  baseline  case,  a  round 
corner  channel  with  1  mm  x  1  mm  square  cross  section  and  4  mm 
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in  length  and  the  rib  with  1  mm  in  width  are  applied.  Five  water 
inlet  pores  with  0.2  mm  diameter  are  attached  to  the  bottom  wall 
of  the  channel  (the  interface  between  a  channel  and  a  GDL).  The 
length  of  the  channel  and  the  size  and  the  location  of  the  pores  are 
determined  according  to  experimental  observations  [46,50-52].  A 
computational  mesh  consisting  of  21 ,460  cells  is  used  for  the  base¬ 
line  case.  The  adequacy  of  this  grid  was  tested  by  increasing  and 


decreasing  the  number  of  grid  nodes  by  ±20%,  and  similar  transport 
processes  were  obtained. 

2.3.  Boundary  and  initial  conditions 

In  the  baseline  case  simulations,  an  air  velocity  of  10  m  s-1  nor¬ 
mal  to  the  boundary  at  0  bar  (gauge  pressure)  and  a  water  velocity 


16ms 


Fig.  2.  Water  motion  in  the  channel  in  the  baseline  case. 
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20ms 


40ms 

a)  Channel  with  a  round  corner 

Fig.  2.  ( Continued ). 


of  1  m  s-1  normal  to  the  boundary  are  applied.  The  operating  tem¬ 
perature  is  set  to  300  K.  The  static  contact  angle  of  the  bottom  wall 
of  the  channel  is  set  to  1 35°,  standing  for  the  typical  PTFE  treatment 
of  a  GDL.  Velocities  of  air  and  water  are  also  set  corresponding  to  the 
normal  vehicular  PEMFC  stacks  condition  to  simulate  real  operating 
environments  [45,46]. 

3.  Results  and  discussion 

Water  dynamic  behaviors  in  a  cathode  serpentine  channel 
with  five  water  inlet  pores  are  simulated.  Effects  of  hydrophilic- 
ity/hydrophobicity  of  channel  walls,  corner  shapes  and  velocities 
of  air  and  water  on  the  water  transport  are  discussed.  The  values 
of  the  static  contact  angles  of  the  channel  walls  are  set  to  45°,  90° 
and  135°,  respectively,  in  the  cases. 

3.1.  Baseline  case 

3.1.1.  Water  motion  in  the  channel  with  a  round  corner 

3. 1.1.1.  Static  contact  angle  45° .  Fig.  2  shows  water  dynamic  behav¬ 
iors  in  the  channel  with  hydrophilicity/hydrophobicity  walls  at 
different  times  (From  left  to  right,  the  static  contact  angles  are 
set  to  45°,  90°  and  135°,  respectively).  It  can  be  seen  that  water 


motions  are  similar  at  the  first  stage  in  all  channels  (0-6  ms)  but 
differ  at  8  ms.  At  the  beginning,  water  moves  toward  one  side  of 
the  channel  due  to  the  blockage  effect  of  the  corner.  Like  that  in  a 
straight  channel  [46],  the  water  is  tore  into  droplets  before  it  moves 
to  the  next  water  injection  pore  because  of  the  hydrophobicity  of 
the  bottom  wall  and  the  shearing  force  of  the  inlet  air.  When  these 
newly  formed  water  droplets  move  to  the  next  water  injection  pore, 
they  collide  with  the  water  droplet  emerging  and  are  pushed  to  the 
outside  wall  of  the  channel  by  the  shearing  force  of  the  air.  Coun¬ 
teracted  by  the  collision  force  and  the  air  shearing  force,  the  surface 
tension  of  the  water  cannot  maintain  the  shape  of  the  water  as  a 
droplet  anymore.  Water  droplets  are  dragged  into  long,  thin  water 
films  and  ascend  to  the  boundary  of  the  side  wall  and  the  top  wall 
of  the  channel.  These  long,  thin  water  films  merge  into  long,  thick 
films  and  move  along  the  boundary  at  1 6  ms.  Water  emerging  from 
pores  at  inlet  part  of  the  channel  ascends  to  the  boundary  of  out¬ 
side  wall  and  the  top  wall  and  eventually  moves  along  the  boundary 
because  of  the  inertia  force  and  the  hydrophilicity  of  the  outside  and 
top  wall.  On  the  contrary,  water  emerging  from  the  corner  and  the 
outlet  part  of  channel  ascends  to  the  boundary  of  the  inner  side  wall 
and  the  top  wall  of  the  channel.  Water  emerging  from  the  inlet  part 
of  the  channel  begins  to  drain  out  of  the  channel  at  18  ms.  Water 
motions  in  the  channel  are  almost  the  same  after  16  ms. 
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Fig.  2.  ( Continued ). 


3. 1.1. 2.  Static  contact  angle  90°  and  135°.  Unlike  those  in  the  chan¬ 
nel  with  hydrophilic  walls,  water  shapes  in  the  channel  with  walls 
with  the  static  contact  angel  set  to  90°  are  mainly  in  the  form  of 
droplets.  In  the  channel  with  hydrophilicity  side  and  top  walls, 
contact  surface  between  water  droplets  and  the  side  and  top  walls 
begins  to  expand  after  6  ms  because  of  the  hydrophilicity  of  the 
walls.  However,  in  the  channel  with  walls  with  the  static  contact 


angel  set  to  90°,  the  surface  tension  of  water  maintains  the  shape 
of  droplets  because  the  hydrophilicity  force  of  the  walls  is  less  than 
the  surface  tension  of  water.  When  the  water  droplet  moves  to  the 
next  water  inject  pore,  it  collides  with  the  water  emerging  from  the 
pore  and  become  a  bigger  one.  Bigger  water  droplets  begin  to  form 
after  16  ms.  Adhering  to  the  outside  wall  by  the  inertial  force  and 
the  air  shearing  force,  water  droplets  emerging  from  the  inlet  part 
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40ms 

b)  Channel  with  a  sharp  corner 


Fig.  2.  ( Continued ). 


of  the  channel  drain  out  of  the  channel.  Although  water  droplets 
emerging  in  the  corner  and  the  outlet  part  of  channel  adhere  to 
the  inside  wall  of  channel,  dynamic  behaviors  of  water  droplets 
emerging  both  in  the  corner  and  outlet  part  and  in  the  inlet  part 
of  the  channel  are  almost  the  same.  Water  motion  in  the  channel 
with  walls  with  the  static  contact  angel  set  to  135°  is  similar  to 


that  in  the  channel  with  walls  with  the  static  contact  angel  set  to 
90°. 

3.1.2.  Water  motion  in  the  channel  with  a  sharp  corner 

3. 1.2.1.  Static  contact  angle  45°.  Duo  to  the  blockage  of  the  sharp 

corner,  the  pressure  in  the  corner  part  of  the  channel  with  a  sharp 
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corner  is  greater  than  that  in  the  channel  with  a  round  corner.  At  the 
very  beginning,  the  water  behavior  in  the  channel  with  a  sharp  cor¬ 
ner  is  similar  to  that  in  the  channel  with  a  round  corner,  as  shown 
in  Fig.  3.  Water  emerging  in  the  inlet  part  of  the  channel  moves  to 
the  outside  wall  of  the  channel  and  finally  forms  a  long,  thin  water 
film  under  the  effect  of  the  hydrophilicity  of  the  outside  wall  and 
the  shearing  force  of  the  inlet  air  at  8  ms.  With  the  further  water 
generation,  water  films  moving  along  outside  wall  ascend  to  the 
boundary  of  the  outside  wall  and  the  top  wall  under  the  shearing 
force  of  the  inlet  air.  When  the  water  films  move  to  the  corner, 
water  emerging  at  the  corner  part  merges  with  these  water  films, 
so  does  the  water  emerging  in  the  outlet  part  of  the  channel.  Since 
the  water  film  is  long  and  relatively  thin,  the  hydrophilicity  of  the 
top  wall  and  the  air  shearing  force  are  greater  than  the  gravity  of 
the  water  film.  So  the  long,  thin  water  films  can  move  on  the  top 
wall  of  the  channel.  Finally,  water  films  move  along  and  cover  the 
whole  top  wall  from  the  corner  part  to  the  outlet  part  of  the  channel 
at  14  ms. 

3.1. 2.2.  Static  contact  angle  90°.  With  the  decrease  of  hydrophilic¬ 
ity  of  walls,  water  emerging  from  the  pores  forms  water  droplets 
rather  than  water  films.  When  water  droplets  moves  to  the  sharp 
corner,  the  blockage  of  the  sharp  corner  and  these  water  droplets 
themselves  forms  a  high  pressure  area  in  the  corner  part.  Water 
droplets  in  the  inlet  part  of  the  channel  move  slowly  and  finally 
merge  into  one  long,  thick  water  film  at  12  ms.  When  the  water 
film  moves  through  the  corner  into  the  outlet  part  of  the  channel, 
it  collides  with  water  emerging  in  the  corner  and  the  outlet  part  of 
the  channel  and  forms  a  thicker  film.  Most  space  inside  the  chan¬ 
nel  is  occupied  by  this  long,  thick  water  film.  Water  behavior  in  the 
outlet  part  of  the  channel  is  quite  different  from  that  in  the  inlet 
part.  Water  exists  almost  in  the  form  of  droplets  during  the  whole 
drainage  process. 

3. 1.2.3.  Static  contact  angle  135°.  The  water  motion  in  the  channel 
with  walls  with  the  static  contact  angel  set  to  135°  is  similar  to  that 
in  the  channel  with  walls  with  the  static  contact  angel  set  to  90°, 
but  water  droplets  are  bigger  and  water  films  are  thicker.  Long  thick 
water  films  occupy  almost  the  whole  inlet  part  and  the  corner  part 
of  the  channel  while  the  outlet  part  of  the  channel  is  occupied  by 
big  water  droplets.  The  blockage  occurs  more  frequently  because 
of  thick  water  films  in  the  inlet  part  of  the  channel  and  big  water 
droplets  in  the  outlet  part  of  the  channel. 

Fig.  3  shows  the  facet  fraction  of  water  on  the  bottom  wall  of 
channels  with  other  walls  set  to  different  static  contact  angles  and 


with  different  corner  shapers  (‘s’  stands  for  a  sharp  corner  and  ‘r’ 
stands  for  a  round  corner).  It  can  be  seen  that  in  all  cases,  the  facet 
fraction  of  water  on  the  bottom  wall  is  less  in  the  channel  with  a 
round  corner  than  in  the  channel  with  a  sharp  corner.  So  the  reac¬ 
tant  gas  diffuses  better  in  the  round  corner  channel  than  in  the 
sharp  corner  channel  and  water  emerges  more  easily  in  the  round 
corner  channel  than  in  the  sharp  corner  channel.  The  volume  frac¬ 
tion  of  water  inside  the  channel  with  a  sharp/round  corner  can 
be  found  in  Fig.  4.  The  round  corner  channel  also  has  less  water 
volume  fraction  than  does  the  sharp  corner  channel.  Due  to  the 
blockage  of  the  sharp  corner,  water  in  the  inlet  part  and  the  cor¬ 
ner  part  moves  slowly.  So  the  water  droplets  in  the  inlet  part  of  the 
channel  tend  to  join  each  other  to  form  even  bigger  ones.  Due  to  the 
weaker  air  shearing  force,  the  consecutive  water  injection  and  the 
slow  movement  of  water  droplets,  water  droplets  on  the  hydropho¬ 
bic  bottom  wall  in  the  inlet  part  of  the  channel  can  connect  with 
each  other  and  eventually  form  long  thick  water  films,  so  do  the 
water  droplets  in  the  sharp  corner  part.  Unlike  the  sharp  corner,  the 
round  corner  provides  a  good  passage  for  both  the  air  and  the  water. 
Due  the  stronger  air  shearing  force,  the  quick  water  droplets  move¬ 
ment  and  the  hydrophobicity  of  the  bottom  wall,  water  droplets  on 
the  hydrophobic  bottom  wall  are  tore  into  smaller  ones  and  move 
quickly  out  of  the  channel.  It  is  more  obvious  in  the  case  with  the 
air  inlet  velocity  set  to  1 5  m  s-1 ,  as  shown  in  Fig.  5  .  It  can  be  seen 
that  lots  of  smaller  water  droplets  are  formed  by  the  shearing  force 
even  in  the  hydrophilic  channel.  Rather  than  forming  long,  thick 
water  films,  water  films  in  the  hydrophilic  channel  that  attach  the 
wall  by  the  surface  tension  are  cut  by  the  strong  shearing  force  into 
short,  thin  films  and  some  small  water  droplets  since  the  surface 
tension  can  only  maintain  a  thin  water  film  under  the  increased  air 
shearing  force. 

However,  it  is  interesting  to  found  that  the  volume  fraction  of 
water  inside  the  channel  decreases  with  the  increase  of  the  static 
contact  angle  in  the  cases  studied  but  the  facet  fraction  of  water 
on  the  bottom  wall  of  the  channel  increase  with  the  increase  of 
the  static  contact  angle.  This  finding  suggests  that  unlike  water  in 
the  channel  with  hydrophobic  walls,  water  in  the  channel  with 
hydrophilic  walls  locates  inside  the  chamber  of  the  channel,  on 
the  top  wall  of  the  channel  or  on  the  side  walls  of  the  chan¬ 
nel  rather  than  on  the  bottom  wall  of  the  channel.  Although  the 
channel  with  hydrophobic  walls  has  the  lowest  volume  fraction  of 
water,  more  water  droplets/films  move  on  the  bottom  wall  than  do 
water  droplets/films  in  the  channel  with  hydrophilic  walls.  With 
the  increase  of  the  air  and  water  inlet  velocity,  the  difference  of 
volume  fraction  of  water  inside  the  channel  of  the  case  studied 
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Fig.  5.  Water  motion  in  the  channel  under  different  air  inlet  velocities. 
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decreases  while  the  difference  of  facet  fraction  of  water  on  the  bot¬ 
tom  wall  maintains.  So  it  may  be  helpful  for  fuel  cells  operating 
at  high  loads  to  expel  water  if  the  channel  walls  are  hydrophilic 
since  the  more  bottom  wall  areas  are  uncovered  by  water,  the  more 
reactant  gas  can  pass  from  the  GDL  to  the  channel. 

3.2.  Effect  of  the  air  inlet  velocity 

Since  the  round  corner  channel  performs  better  than  the  sharp 
corner  channel,  the  discussion  below  is  based  on  the  round  corner 
channel  model. 

3.2.1.  Static  contact  angle  45° 

The  air  inlet  velocity  has  a  significant  effect  on  the  dynamic 
behavior  of  water.  Simulation  cases  in  which  the  air  inlet  velocity 
set  to  5, 10  and  15  m  s-1  are  studied  to  investigate  water  transport 
in  the  channel. 

Results  show  that  water  motions  at  different  air  inlet  velocities 
are  quite  different.  Fig.  5  shows  water  motions  at  different  air  inlet 
velocities.  When  the  inlet  air  velocity  equals  5  m  s-1 ,  water  droplets 
emerging  from  inlet  pores  keep  the  shape  of  droplets  because  the 
shearing  force  of  air  is  less  than  the  surface  tension  of  the  water. 

Water  droplets  are  elongated  at  a  larger  air  inlet  velocity.  Water 
droplets  attach  to  the  sidewall  of  the  channel  at  4  ms  because  the 
hydrophilicity  of  the  wall  is  greater  than  the  shearing  force  of  the 
air.  Water  emerging  from  the  inlet  part  of  the  channel  attaches  to 
the  outside  wall  at  the  very  beginning  and  moves  to  the  boundary 
part  of  the  outside  wall  and  the  top  wall.  So  does  water  emerging  % 

from  the  corner  part.  Finally,  water  emerging  from  the  inlet  part  £ 

and  the  corner  of  the  channel  merges  into  one  long,  thin  water  film 
and  moves  along  the  boundary  of  the  outside  wall  and  the  top  wall. 

On  the  contrary,  water  emerging  from  the  outlet  part  of  the  channel 
moves  towards  the  boundary  of  the  inside  wall  and  the  top  wall  and 
then  drain  out. 

When  the  inlet  air  velocity  reaches  1 5  m  s-1 ,  the  shearing  force 
of  the  air  is  greater  than  that  in  the  baseline  case,  so  water  emerg¬ 
ing  at  the  inlet  part  is  tore  into  small  water  droplets.  Small  water 
droplets  collide  with  each  other  and  form  several  thin  water  films 
under  the  air  shearing  force  and  hydrophilicity  of  channel  walls. 
Dragged  by  the  inertial  force,  those  water  films  form  one  long, 
thin  film  when  they  move  through  the  corner.  Water  emerging  in  | 
the  corner  is  tore  into  many  tiny  water  droplets  and  are  blew  out  [o 

through  the  chamber  without  attaching  to  any  wall  of  the  chan-  | 

nel.  Unlike  that  of  water  emerging  in  the  inlet  part  and  the  corner,  ^ 
motion  of  water  emerging  in  the  outlet  part  of  the  channel  is  mainly 
in  the  form  of  thick  water  films. 

3.2.2.  Static  contact  angle  90° 

Compared  to  that  in  the  baseline  condition,  the  air  shearing 
force  is  less  when  the  inlet  air  velocity  equals  5 ms-1.  Under  a 
less  shearing  force,  water  emerging  from  the  pores  is  inclined  to 
form  bigger  water  droplets  under  the  surface  tension.  It  can  be  seen 
that  when  water  droplets  move  through  the  water  inlet  pores,  they 
merge  with  the  water  emerging  from  the  pore  and  form  even  big¬ 
ger  ones.  Water  motion  in  the  channel  is  mainly  in  the  form  of  big 
water  droplets.  Less  hydrophilicity  and  greater  shearing  force  lead 
to  small  water  droplets  in  the  channel  when  the  air  inlet  veloc¬ 
ity  reaches  1 5  m  s-1 .  Even  though  water  droplets  collide  with  each 
other,  newly  formed  water  droplets  are  so  small  that  they  drain  out 
of  the  channel  quickly. 

3.2.3.  Static  contact  angle  135° 

The  water  motion  is  similar  to  that  in  the  channel  with  the  static 
contact  angle  of  walls  set  to  90°  under  both  conditions.  It  can  be 
seen  that  the  difference  of  water  motions  in  above  cases  diminishes 
with  the  increase  of  the  air  inlet  velocity. 


■  contact  angle  45°(bottom  wall) 
0.5  -i  •  contact  angle  45°(top  wall) 
-A-  contact  angle  45°(volume) 


— □ —  contact  angle  90°(bottom  wall) 
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Fig.  6.  Facet  and  volume  fraction  of  water  in  the  channel  under  different  air  inlet 
velocities. 

Fig.  6  shows  the  volume  fraction  of  water  in  the  channel  and  the 
facet  fraction  of  water  on  the  walls  of  the  channel  under  different 
air  velocities.  It  can  be  seen  that  the  volume  fraction  of  water  inside 
the  channel  drops  with  the  increase  of  the  static  contact  angle.  The 
drop  is  noticeable  when  the  air  velocity  increases  from  5  m  s-1  to 
10ms_1  but  less  significant  when  the  air  velocity  increases  from 
10  ms-1  to  15  ms-1.  It  also  can  be  found  that  under  different  air 
velocities,  the  channel  with  wall  static  contact  angle  set  to  45°  has 
the  least  facet  fraction  of  water  on  the  bottom  wall  but  the  high¬ 
est  volume  fraction  of  water  inside  the  channel.  Since  the  shearing 
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Fig.  7.  Water  motion  in  the  channel  under  different  water  inlet  velocities 
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force  of  the  air  is  not  great  enough  to  counteract  the  effect  of  the 
surface  tension  of  the  water,  water  merges  into  thick  films  or  big 
droplets  under  air  velocity  at  5  m  s-1 .  Water  films  and  droplets  are 
tore  apart  with  the  increase  of  the  inlet  air  velocity.  These  thin, 
short  films  and  small  droplets  are  easier  to  drain  out  than  are  thick 
films  and  big  droplets. 

3.3.  Effect  of  the  water  inlet  velocity 

3.3. 1.  Static  contact  angle  45° 

Fig.  7  shows  the  water  motion  under  different  water  inlet  veloc¬ 
ities.  Under  less  water  generating  rate  (0.5  m  s-1 ),  water  emerging 
from  the  pores  is  tore  into  many  small  water  droplets.  Unlike  those 
in  the  baseline  case,  water  droplets  under  this  circumstance  collide 
and  form  several  water  films  rather  than  one  long  water  film.  Water 
emerging  from  the  outlet  of  the  channel  exists  mainly  in  the  form 
of  small  water  droplets.  Water  forms  small  water  droplets  at  a  high 
water  generating  rate  (1.5  m  s-1 ).  Water  droplets  collide  and  form 
many  big  water  droplets  rather  than  one  water  film. 

3.3.2.  Static  contact  angle  90° 

The  water  motion  in  the  channel  with  water  generating  rate 
at  0.5 ms-1  is  also  different  from  that  in  the  baseline  case.  Water 


forms  droplets  more  quickly  than  that  in  the  baseline  case.  The  size 
of  the  water  droplets  formed  under  collision  in  the  inlet  part  of  the 
channel  is  small.  These  water  droplets  move  through  the  corner  on 
the  outside  wall  of  the  channel  and  finally  drain  out.  When  water 
generating  rate  equals  1.5 ms-1,  the  water  motion  in  the  channel 
with  the  wall  static  contact  angle  set  to  90°  is  almost  the  same  as 
that  in  the  channel  with  the  wall  static  contact  angle  set  to  45°. 

3.3.3.  Static  contact  angle  135° 

Under  this  condition,  water  in  the  channel  forms  more  water 
droplets  at  an  earlier  stage.  Water  drains  out  mainly  in  the  form  of 
small  droplets.  When  water  generating  rate  equals  1.5 ms-1,  the 
water  motion  in  the  channel  with  the  wall  static  contact  angle  set 
to  1 35°  is  almost  the  same  as  that  in  the  channel  with  the  wall  static 
contact  angle  set  to  45°  and  90°  except  for  the  size  of  the  droplets. 
Because  of  the  hydrophobicity  of  walls,  water  mainly  moves  on  the 
bottom  wall  rather  than  on  the  side  walls. 

Fig.  8  shows  the  facet  and  volume  fraction  of  water  in  the  chan¬ 
nel  under  different  water  inlet  velocities.  It  is  interesting  to  found 
that  under  all  the  conditions,  the  channel  with  water  inlet  veloc¬ 
ity  at  1  ms-1  has  the  highest  water  fraction  on  the  bottom  wall 
while  the  channel  with  the  least  water  inlet  velocity  has  the  least 
water  volume  fraction.  When  the  water  inlet  velocity  is  low,  water 


Facet/Volume  fraction  Facet/Volume  fraction  Facet/Volume  fraction 
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Fig.  8.  Facet  and  volume  fraction  of  water  in  the  channel  under  different  water  inlet 
velocities. 


entering  the  channel  forms  only  small  droplets  or  thin  water  films 
and  the  inertia  force  is  not  strong  enough  to  push  these  droplets 
or  films  to  the  side  wall  or  top  wall.  So  the  water  droplets  and 
films  move  through  the  channel  mainly  on  the  bottom  wall.  The 
volume  fraction  of  water  is  less  because  water  entering  the  chan¬ 
nel  in  certain  time  is  less  due  to  the  low  water  inlet  velocity  and 
because  small  droplets  and  thin  films  are  easy  to  drain  out  of  the 
channel.  On  the  contrary,  water  forms  larger  droplets  or  thicker 
films  when  water  inlet  velocity  increases.  Under  this  circumstance, 
large  water  droplets  or  thick  water  films  have  more  kinetic  energy 
and  are  pushed  toward  the  side  wall  or  even  the  top  wall.  So  water 
moves  mainly  on  the  side  wall  or  on  the  top  wall.  Since,  in  a  cer¬ 
tain  time,  more  water  enters  the  channel  under  a  high  water  inlet 
velocity  and  large  droplets  and  thick  films  are  hard  to  drain  out  of 
the  channel,  the  volume  fraction  in  the  channel  with  a  high  water 
inlet  velocity  is  greater  than  that  in  the  channel  with  a  low  water 
inlet  velocity. 

Besides,  the  pressure  loss  in  the  channel  is  investigated  and 
findings  show  that  pressure  loss  increases  with  the  increase  of 
hydrophobicity  of  channel  walls.  In  channel  with  hydrophobic- 
ity  walls,  water  occurs  mainly  in  the  form  of  droplets  as  seen  in 
Figs.  2,  5  and  7.  Unlike  the  water  films  shaped  in  the  channel  with 
hydrophilic  walls,  water  droplets  in  the  channel  with  hydrophobic 
walls  need  more  energy  to  move,  resulting  in  more  pressure  loss.  So 
water  droplets  in  the  channel  with  hydrophobic  walls  move  out  of 
the  channel  more  quickly  than  the  water  films  and  water  droplets  in 
the  channel  with  hydrophilic  walls  in  the  cost  of  more  pressure  loss. 
The  results  match  those  in  both  experimental  and  computational 
works  [12,13,43]. 

4.  Conclusions 

The  following  conclusions  can  be  obtained  from  above  study. 

1 )  Water  drains  out  more  quickly  in  a  round  corner  channel  than  in 
a  sharp  corner  channel.  Water  moves  in  the  shape  of  thick  films 
in  the  inlet  and  corner  part  of  a  sharp  corner  channel  but  in  the 
shape  of  droplets  in  the  inlet  and  corner  part  of  a  round  corner 
channel. 

2)  The  volume  fraction  of  water  inside  the  channel  decreases  with 
the  increase  of  the  static  contact  angle  but  the  facet  fraction 
of  water  on  the  bottom  wall  of  the  channel  increase  with  the 
increase  of  the  static  contact  angle.  The  formation  and  motion 
of  water  droplets/films  in  a  channel  with  a  corner  are  quite  dif¬ 
ferent  from  those  in  a  straight  channel. 

3)  Water  emerging  from  one  pore  has  a  great  impact  on  the  for¬ 
mation  of  and  subsequent  on  the  behavior  of  water  emerging 
from  other  pores.  The  formations  and  movements  of  water  dif¬ 
fer  under  different  air  inlet  velocities,  water  inlet  velocities  and 
wall  static  contact  angles,  but  the  difference  diminishes  with  the 
increase  of  the  air  inlet  velocity. 
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